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Abstract: 

Fuzzy spheres appear as classical solutions in a matrix model obtained via dimensional 
reduction of 3-dimensional Yang-Mills theory with the Chern-Simons term. Well-defined 
perturbative expansion around these solutions can be formulated even for finite matrix 
size, and in the case of k coincident fuzzy spheres it gives rise to a regularized U(A;) gauge 
theory on a noncommutative geometry. Here we study the matrix model nonperturbatively 
by Monte Carlo simulation. The system undergoes a first order phase transition as we 
change the coefficient (a) of the Chern-Simons term. In the small a phase, the large 
properties of the system are qualitatively the same as in the pure Yang-Mills model (a = 0), 
whereas in the large a phase a single fuzzy sphere emerges dynamically. Various 'multi 
fuzzy spheres' are observed as meta-stable states, and we argue in particular that the k 
coincident fuzzy spheres cannot be realized as the true vacuum in this model even in the 
large N limit. We also perform one-loop calculations of various observables for arbitrary 
k including k = \. Comparison with our Monte Carlo data suggests that higher order 
corrections are suppressed in the large N limit. 
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1. Introduction 

Quantization of gravity is one of the most important problems in particle physics, and lots 
of attempts have been made so far. String theory, in particular, enables us to quantize small 
fluctuations of the space-time metric around certain backgrounds, but the nonperturbative 
determination of the background itself is yet to be done. Such an issue may be addressed 
by studying matrix models, which are proposed as nonperturbative formulations of string 
theory. For instance, the IKKT matrix model [1,2], which can be obtained via dimensional 
reduction of 10-dimensional super Yang-Mills theory, is conjectured to be a nonperturbative 
definition of type IIB superstring theory in ten dimensions. 

In this model the space-time is represented by the eigenvalue distribution of ten bosonic 
matrices, and if the distribution collapses to a four-dimensional hypersurface, we may 
naturally understand the dimensionality of our space-time as a result of the nonperturbative 
dynamics of superstring theory [3]. Prom the path-integral point of view, this phenomenon 
may be caused by the phase of the fermion determinant [4] , and Monte Carlo results support 
this mechanism [5] . In Ref . [6] the first evidence for the above scenario has been obtained by 
calculating the free energy of the space-time with various dimensionality using the Gaussian 
expansion method up to the 3rd order. Higher-order calculations [7] and the tests of the 
method itself in simpler models [8] have strengthened the conclusion considerably. Another 
indication of this phenomenon has been obtained recently from the calculations of the 2- 
loop effective action around fuzzy spheres [9]. See also Refs. [10-14] for related works. 

Since the space-time is described by matrices in the IKKT model, it is generically 
noncommutative. If one expands the model around a D-brane background for instance, 
one obtains a gauge theory on noncommutative (NC) geometry [15]. ^ This connection 
was also understood solely in terms of string theory [19], which triggered much interest in 
the dynamical properties of field theories on NC geometry in general. It was found that 
various interesting phenomena which do not have the commutative counterparts occur due 
to the so-called UV/IR mixing effects [20]. By considering finite-dimensional matrices, one 
may naturally regularize NC field theories. In particular the twisted reduced models [21], 
which appeared in history as an equivalent description [22] of large N gauge theories, can 
be interpreted as a lattice regularization of field theories on a noncommutative torus [23]. 
This enables nonperturbative studies of various dynamical issues in these theories by Monte 
Carlo simulation [24-26]. See Ref. [27] for a general review on the dynamics of matrix 
models in the context of superstring theory and noncommutative geometry. 

A different, although closely related, type of regularization is known as the 'fuzzy 
sphere' [28], which can also be formulated using finitc-dimcnsional matrices. The UV regu- 
larization in this case is introduced by putting an upper bound on the angular momentum 
when one expands a function on the sphere in terms of spherical harmonics. This cutoff 
procedure is compatible with the 'star-product' of functions which appears in NC geometry. 
Since the regularized theory does not break the continuous symmetries of the sphere, the 

^For an earlier work on the connection between matrix models and NC geometry in the context of 

toroidal compactification, see [16]. A mathematical formulation of field theories on NC geometry is given 
in Ref. [17]. For a comprehensive review including recent developments, see Ref. [18]. 
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well-known problems in lattice field theory concerning chiral fermions and super symmetry 
may be easier to overcome [29]. The first challenge in this direction is to remove the UV/IR 
mixing effects in the 'continuum limit'. 

Being the simplest curved compact noncommutative space, the fuzzy 2-sphere has been 
studied extensively in the literature. One of the fundamental issues is the construction of 
the Dirac operator. The one proposed in Refs. [30] does not have the fermion doubling 
problem [31,32] and it yields the correct chiral anomaly both in the global form [33-35] and 
in the local form [36]. This Dirac operator, however, breaks chiral symmetry explicitly, 
and it may therefore be regarded as an analog of the Wilson fermion in lattice gauge 
theory. On the other hand, an analog of the overlap Dirac operator [37], which satisfies 
the Ginsparg- Wilson relation [38] and hence preserves modified chiral symmetry [39], is 
constructed for the free fermion [40] and for general gauge configurations [41,42]. ^ In 
this case the chiral anomaly arises from the measure and the correct results are reproduced 
in Refs. [41,48,49]. Gauge configurations with non-trivial topology have been found for 
the former Dirac operator [34,50-52] as well as for the latter [48,53]. The Seiberg-Witten 
map [19] has been constructed also on the fuzzy 2-sphere [54]. 

In string theory fuzzy spheres appear as classical solutions in the presence of an external 
Ramond-Ramond field [55], and the low-energy effective theory is given by NC Yang-Mills 
theory with the Chcrn-Simons term [56] (Sec also [57]). Classical solutions of the effective 
theory and their D-brane interpretation have been studied in Ref. [58]. In the matrix-model 
description of string theory, fuzzy spheres appear as solutions to the classical equation 
of motion if one adds the Chern-Simons term representing the coupling to the external 
field [59] . Expanding the bosonic matrices around the classical solution, one obtains a NC 
gauge theory on fuzzy spheres. The situation is quite analogous to the flat D-brane in 
matrix models [15] except that the fuzzy sphere can be realized even for finite matrices, 
which makes various calculations totally well-defined. Studying fuzzy spheres is expected 
to give us a new insight into the description of curved space-time in matrix models, which is 
vitally important in the context of quantum gravity. A variety of fuzzy-sphere-like solutions 
in matrix models have been studied in Refs. [9,59-65]. Fuzzy spheres are also discussed in 
a matrix model [66] for the M-theory in the so-called pp-wave background [67]. 

The stability of the fuzzy sphere against quantum fluctuations is a non-trivial impor- 
tant issue, which has been discussed by perturbative calculations [9,59,64,68] and by the 
Gaussian expansion [65]. Since matrix models typically have various fuzzy-sphere-like so- 
lutions, it is important to determine which one describes the true vacuum. Such an issue 
is related to the dynamical generation of not only the space-time but also the gauge group, 
since k coincident fuzzy spheres give rise to a NC gauge theory with the gauge group of 
rank k. The unifled treatment of the space-time and the gauge group is one of the advan- 
tages of considering noncommutative geometry or the matrix model formulation of string 

^The overlap Dirac operator for general gauge backgrounds has been proposed earlier on a NC torus 
[43], and used for a construction of chiral gauge theories. The correct chiral anomaly in this case has 
been reproduced in Ref. [44] for arbitrary even dimensions by using a topological argument [45]. In odd 
dimensions the analogous Dirac operator is used to define a lattice NC Chern-Simons theory [46] through 
the parity anomaly [47]. 
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theory. 

In this paper we study a matrix model obtained via dimensional reduction of three- 
dimensional Yang-Mills-Chern-Simons theory, which is known to have fuzzy spheres as 
classical solutions. Unlike previous works we perform fully nonperturbative first-principle 
calculations using Monte Carlo simulations. We find that the single fuzzy sphere is non- 
perturbatively stable if the coefficient (a) of the Chern-Simons term is sufficiently large, 
but it collapses to a 'solid ball' if a is smaller than a critical value. The transition be- 
tween the two phases is of first order, and we observe a strong hysteresis. Various 'multi 
fuzzy spheres', which are also classical solutions of the model, are observed as meta-stable 
states, and we argue in particular that the k coincident fuzzy spheres cannot be realized 
as the true vacuum even in the large A'^ limit. We also perform one-loop calculations of 
various observables for arbitrary k including k = I. Comparison with our Monte Carlo 
data suggests that higher order corrections are suppressed in the large N limit. 

The rest of the paper is organized as follows. In Section 2 we define the model and 
review how fuzzy spheres appear in this model as classical solutions. In Section 3 we 
investigate the phase diagram of the model and demonstrate the existence of the first 
order phase transition. In Section 4 wc discuss the properties of the single fuzzy sphere. In 
Section 5 we study the geometrical structure of the dominant configurations in each phase. 
In Section 6 we show how various multi fuzzy spheres appear as meta-stable states and 
study their properties. Section 7 is devoted to a summary and discussions. In Appendix 
A we comment on the algorithm used for our Monte Carlo simulations. In Appendices B 
and C we give a self-contained derivation of the one-loop results for the effective potential 
as well as various observables. In Appendix D we discuss the instability of the k coincident 
fuzzy spheres. 



2. The Yang-Mills-Chern-Simons matrix model and the fuzzy spheres 
The model we study in this paper is defined by the partition function [59] 

Z = j dAfT^ , (2.1) 

S = Ntx (^-^ [Ai„A^]^ + ^iaei,^pA^A^Ap^ , (2.2) 

where A^ {fi = 1,2,3) are N x N traceless Hermitian matrices. The integration measure 
dA is defined by d^ = Ua=i U^i ^> where Af^ are the coefficients in the expansion 
Af^ = X^a^Jii" with respect to the SU(A^) generators normalized as tr{t"-t^) = ^6°''^. 
The coefficient of the first term in (2.2) is fixed, but this does not spoil any generality since 
one can always rescale A^ to bring the action into the present form. 

This model may be regarded as the zero- volume limit of SU{N) Yang- Mills theory with 
the Chern-Simons term in the three-dimensional Euclidean space, and it has the SO (3) 
rotational symmetry as well as the SU(A^) symmetry A^ U A^U'^ , where U G SU(A^). 
There are some points to note here, though. The parameter a in (2.2) is chosen to be real in 
order for fuzzy spheres to be classical solutions of the model. As a result the Chern-Simons 



-4- 



term in the action (2.2) is real unlike the counterpart in ordinary field theory, where it is 
purely imaginary and therefore poses a severe technical problem in Monte Carlo simulation. 
Note also that a may take arbitrary real number without breaking any symmetries unlike 
the ordinary Chern-Simons theory, where the coefficient is quantized for the invariance 
under topologically nontrivial gauge transformations. Since the model has the duality 
OL 1-^ —OL associated with the parity transformation i— —A^)^^ we restrict ourselves to 
a > throughout this paper without loss of generality. 

The pure Yang-Mills model {a = 0) and its obvious generalization to D dimensions 
with D matrices (/x = 1, • • • , have been studied by many authors. In particular the 
large N dynamics of the model have been studied by the 1/D expansion and Monte Carlo 
simulations [69]. The partition function was conjectured [71] and proved [72] to be finite for 
N > D /{D — 2). (Sec Rcfs. [73, 74] for the supersymmetric case.) The partition function in 
the presence of the Chern-Simons term has been studied analytically for N = 2 [75], and 
it turned out to be convergent in the supersymmetric case, but not in the bosonic case. It 
is also proved that adding a Myers term (the Chern-Simons term in the present case) does 
not affect the convergence as far as the original path integral converges absolutely [76], 
which means in particular that the partition function (2.1) is convergent for iV > 4. Note 
that this statement holds despite the fact that the classical action (2.2) is unbounded from 
below for a 7^ 0. 

The classical equation of motion in the present model reads 

[A,,Ap]=0. (2.3) 

The simplest type of solutions is given by the commutative matrices, which can be brought 
into the diagonal form 

A^ = diag(x«,42),... ,a;(f)) (2.4) 

by an appropriate SU(A^) transformation. This type of solutions exists also at a = 0, and 
the one-loop efi'ective action around it has been calculated in Ref. [69]. For a 7^ the 
one-loop effective action reads [59] 



Tl-loop — ^ 

i<j 



log|(x»-x(f))'Ulog|l-— 



(2.5) 



where the second term represents the additional piece for a 7^ 0. The resulting attrac- 

(i) 

five force between xli makes their distribution shrink until the perturbative calculation 
becomes no more valid [69]. 

When a is nonzero, there exists another type of solutions given by 

A^ = aLi,, (2.6) 

where (/x = 1, 2, 3) satisfy 

[Ln,L^] = iCfj^ypLp . (2-7) 
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Since the equation (2.7) is nothing but the SU(2) Lie algebra, each of its A^-dimensional 
representation yields a classical solution. We denote the n-dimensional irreducible repre- 
sentation of the algebra by L^T^ , which can be characterized by the identity 



E(4"'r = i("^-i)i»- 



(2.8) 



Let us consider the case in which itself is given by the irreducible representation. 
The corresponding classical solution 



satisfies 



R^l 



N 



(2.9) 



(2.10) 



R 



■ a 



(2.11) 



If we neglect for the moment the fact that A^^ are noncommutative, eq. (2.10) implies that 
all the N eigenvalues of Af^ arc distributed on a sphere centered at the origin with the 
radius R. Since A^ are noncommutative in reality, this classical solution is called the fuzzy 
2-sphere (In this paper wc refer to it as the ^single fuzzy sphere'). 

If we consider a general representation, the corresponding classical solution can be 
brought into the form 



A,^ 



a 



L 



("2) 



V 



by an appropriate SU(A'^) transformation, where 

k 



a=l 



N 



This solution satisfies 



/( 



n )2 1 



{r2? In. 



where 



1 



\ 



r(™fc) / 



1 



(2.12) 



(2.13) 



(2.14) 



(2.15) 



meaning that Ua eigenvalues of Afj_ are distributed on a sphere centered at the origin with 
the radius Tq. In the k > 2 case we call the solution (2.12) ^multi fuzzy spheres'. 
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The particular case of 'multi fuzzy spheres' that we will focus on later corresponds to 
taking rai = • • • = nfc = n, and therefore N = n-k due to (2.13). The radius (2.15) for each 
sphere becomes equal, and it is proportional to 1/k at large N. Expansion of the model 
around such a configuration, which describes k coincident fuzzy spheres, gives rise to a 
NC gauge theory on the fuzzy sphere with the gauge group of rank k [59] . This statement 
holds also for k = I, which corresponds to the single fuzzy sphere. 

We study the matrix model (2.2) by Monte Carlo simulations using the heat-bath 
algorithm, which was adopted also in the study of the pure Yang-Mills model {a = 0) in 
Ref. [69]. See Appendix A for more details. 



3. The phase structure 

In this Section we study the phase structure of the model by calculating fundamental 
quantities such as (S) and (;^tr(A^)^) by Monte Carlo simulations. Here we take the 
initial configuration to be either of the following two. 



a L 




(N) 



(the single fuzzy sphere start) 
(the zero start) 



(3.1) 



Fig. 1 shows the results plotted against a for iV = 8, 16, 24. The one-loop results obtained 
from the perturbative expansion around the single fuzzy sphere are plotted for comparison. 
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Figure 1: The observables -^{S) (left) and (-^trA^) (right) are plotted against a for A'' = 8, 16, 24. 
The open symbols represent the single fuzzy sphere start, whereas the closed symbols represent the 
zero start. We also plot the one-loop results for the single fuzzy sphere at each N for comparison. 

The Monte Carlo results depend on the initial configuration in the intermediate region 
of a for A?^ > 16, and we observe discontinuities at 

{air ~ (for the single fuzzy sphere start) 

Ocr^ ^ 0.66 (for the zero start) , 

which we call the lower/upper critical points, respectively. This clearly demonstrates the 
existence of a first order phase transition. 
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Note that the large N behavior of the lower critical point is different from the upper 
critical point. If we naively apply the formula (3.2) to = 8, we obtain ac) — ^ — 0.742, 

which is actually larger than Ocr ^ • What happens in reality is that the hysteresis simply 
disappears for = 8 as seen in Fig. 1. 

In the small a phase the results do not depend much on a. In particular the large N 
behavior of each observable is qualitatively the same as in the pure Yang-Mills model [69] 
{a = 0), i.e., -^{S) ~ 0(1) and (itr(^^)2) ~ 0(1). We will call this phase the 'Yang- 
Mills phase'. In the large a phase the results agree very well with the one- loop pcrturbative 
calculations around the single fuzzy sphere, which we elucidate further in Section 4. We 
will call this phase the 'fuzzy-sphere phase'. When we perform simulations with the zero 
start for a > we observe various multi fuzzy spheres as meta-stable states in the 

thermalization process (See Section 6.1). 

4. Properties of the single fuzzy sphere 

In this Section we study various properties of the single fuzzy sphere by Monte Carlo 
simulations using (2.9) as the initial configuration. 

4.1 The low^er critical point and the 'one- loop dominance' 

Here we calculate various observables by Monte Carlo simulations, and compare the results 
with the one-loop pcrturbative calculations. We determine the lower critical point, and 
show that our Monte Carlo data above the critical point agree very well with the one-loop 
results. 

Let us first summarize the results of one-loop calculations. For brevity we introduce 
the notation 



so that the action reads S = N'^{-^tvF'^ + aM). The leading large N behavior of various 
observables is given by 



(4.1) 
(4.2) 




(4.5) 



(4.6) 



(4.3) 



(4.4) 



where we have introduced the rescaled parameter 



a = ay/N . 



(4.7) 
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Figure 2: Various observables are plotted against a = a^/N for = 8, 16, 32 with the single fuzzy 
sphere start. The solid lines represent the classical results, whereas the dashed lines represent the 
results including the one- loop corrections. 



In (4.3) ~ (4.6) the first term represents the classical result, while the second term repre- 
sents the one-loop correction (See Appendix C for derivation). 

In Fig. 2 we plot our Monte Carlo re- 
sults for these quantities against a. All the 
observables exhibit a discontinuity at 



2.16 



a(log(N)/NT + 2.086.. 



2.14 



a 



~ 2.1 . 



(4.8) 2.12 



In terms of the original parameter a, it cor- 
responds to the lower critical point given in 
(3.2). 

The results above the lower critical point 
scale with AT, and moreover they agree very 
well with the perturbativc results including 
the one-loop corrections. For (■^tr(F^jy)^), 
the data agree even with the classical result 
due to the accidental cancellation of the one- 
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2.08 



0.004 



0.008 



0.012 



log{N)/N'' 



Figure 3: The lower critical point is plotted 
for N = 16, 20, 24. The dashed line represents a 
(1)^^^ + a where a = 2.781. 



fit to Qcr 
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loop correction (4.5). These results suggest that the higher order corrections axe suppressed 
at large A'' for fixed a. 

In fact it turns out that the lower critical point (4.8) can be reproduced from the 
effective potential for the scalar mode on the fuzzy sphere [70] as ^ 



3/4 



+ a 



logN 



(4.9) 



In Fig. 3 we plot the lower critical point ctcr determined by simulations for N = 8, 16, 24, 
where the 'error bars' represent the range of a for which we observe two-state signals. Our 
results are consistent with the large N asymptotic behavior (4.9). 

4.2 Reproducing an exact result 

In this Section we consider an exact result, which can be used as a check of Monte Carlo 
simulations and the 'one-loop dominance' discussed in the previous Section. Let us note 
that the expectation value of the observable 



K 



+ 3aM 



is given exactly by 



(4.10) 



(4.11) 



for arbitrary a and N. This identity has been studied also in the a = case [69]. 

The exact result (4.11) can be derived 
from the invariancc of the partition func- 
tion (2.1) under the change of integration 



variables Ar 



[1 + e)A^,. Alternatively 



one may derive it from 

d 



o 
V 



0= / d^^{tr(tM^)e-^} ,(4.12) 
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2 
a 



dA 

which is one of the Schwinger-Dyson equa- 
tions in this model. 

In Fig. 4 we show our results for the 
ratio All the data for various a and 

N are consistent with '1' within error bars, 
which demonstrates the validity of our sim- 
ulation. 

The exact result (4.11) is consistent with our assertion that higher loop corrections are 
suppressed at large in the fuzzy sphere phase. Using (4.5) and (4.6), we obtain 

1 



Figure 4: The ratio is plotted against a for 
N = 8, 16, 32 with the single fuzzy sphere start. 
The horizontal line corresponds to the exact re- 
sult. 



{K) 



1— loop 



-trF^ 



N 



- + 3a 



1— loop 
3 



3q;(M)i_ioop 



a" 1 
D a 



3 , 



(4.13) 



^This was informed to us by Denjoe O'Connor after J.N. gave a seminar on the main results of this 
paper including (4.8) at Dublin Institute for Advanced Studies. 
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which agrees with the exact result (4.11) in the large N limit. We also notice that as a 
becomes large, there is a huge cancellation between the first term (-^tr(F^y)^) and the 
second term 3a{M). Namely the large a behavior of each term is 0{a^), but the leading 
term cancels exactly leaving behind an 0(1) quantity. This is reflected in Fig. 4, where the 
error bars become much larger at d > a^r ■ Such a cancellation is absent in the Yang-Mills 
phase since the second term 3a{M) is close to zero as one can see from Fig. 2. 

5. The eigenvalue distribution of the Casimir operator 

In the previous Section we have seen that our Monte Carlo results in the fuzzy sphere phase 

agree very well with the one-loop calculations around the single fuzzy sphere. Here we show 
more directly that the dominant configurations in this phase indeed have the geometry of 
the 2-sphere. 

For that purpose it is convenient to consider the distribution 

1 ^ 

where \j (j = 1, • • • , A^) are N real positive eigenvalues of the Casimir operator 

Q = • (5-2) 

Note that the set of eigenvalues {Aj} is invariant under both S0(3) and SU(A'^) transfor- 
mations. Within the classical approximation, the distribution f(x) is given by 

f{x) = d{x-R') (5.3) 

for the single fuzzy sphere (2.9), where R is defined in (2.11). Similarly for the multi fuzzy 
spheres (2.12) we get 

/(^) = E (77) '^(^ - ^a) , (5.4) 
a=l 

where Va is defined in (2.15). This will play an important role in identifying various multi 
fuzzy spheres in Section 6. 

Fig. 5 shows the results for a = 0.5 < (left) and a = 1.0,2.0 > (right) at 

A'^ = 16. The vertical lines represent the delta function corresponding to the classical result 
(5.3) for the single fuzzy sphere. At a = 1.0 and 2.0 the measured distribution is close to 
the classical result. Slight deviations such as the smearing and the shift can be understood 
as quantum effects, as we will see in Section 5.1. Thus the dominant configurations are 
not given exactly by the single fuzzy sphere (2.9), but they still preserve the geometrical 
structure of a 2-sphere with a slightly smaller radius and with a finite width. 

The result for a = 0.5, on the other hand, deviates drastically from the classical result 
(5.3). In Section 5.2 we will discuss how the observed distribution can be understood 
qualitatively by considering the diagonal configuration (2.4) and the effective action (2.5). 
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Figure 5: The eigenvalue distribution f{x) of the Casimir operator (5.2) is plotted for A'' = 16 at 
a = 0.5 < a^) (left) and a = 1.0, 2.0 > air^ (right). 



At the intermediate values of a we observe a two-state signal, and as a result the mea- 
sured distribution is given by the superposition of the two types of distribution described 
above. This is expected from the first-order nature of the phase transition revealed in 
Section 3. 
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Figure 6: The eigenvalue distribution f{x) of the Casimir operator (5.2) is plotted against x — xo 

for N = 8 (upper left), A'^ = 16 (upper right), A'' = 32 (lower left) and A'' = 64 (lower right), 
respectively. At each A^ the results for different a lie on top of each other. 



- 12 - 



5.1 The distribution in the fuzzy sphere phase 

Here we discuss the properties of the distribution f{x) in the fuzzy sphere phase in more 
detail, and show that they can be understood by the one-loop dominance. 

Let us start with the position of the 
peak, which may be defined by 



Jo 



dxxf{x) = {^triA^f 



(5.5) 




Figure 7: The variance of the distribution f{x) 
is plotted for N = 8, 16, 32, 48, 64. The straight line 
is a fit to (T^ ~ 2 logN + b, where b = -0.84. 



This quantity is actually studied in Sec- 
tion 4.1 and it agrees with the one- loop re- 
sult (4.4) for the fuzzy sphere. Note that 
the one-loop correction is negative, which 
implies that the peak moves toward the 
origin compared with the classical result 

(5.3) . Also the one-loop correction is pro- 
portional to J^, meaning that the shift 
of the peak becomes more pronounced at 
smaller a. These properties are indeed 

observed in Fig. 5 (right). An interesting point to note here is that the one-loop result 

(4.4) for -^(-^tr(^^)) becomes zero at a = ^2. Assuming the one-loop dominance, wc may 
suspect that nonperturbative effects should become non-negligible as we decrease a, which 
invalidates the perturbation theory before d reaches V2. This is indeed what happens in 
reality (Note that a^cr = 2.1 > ^/2). 

Having understood the position of the peak, let us turn to the shape of the distribution. 
In Fig. 6 we plot the distribution f{x) against x — xq for various d at iV = 8, 16,32,64. 
The results for different a collapse to a single curve, which means that the shape of the 
distribution does not depend on a. Furthermore the distribution turns out to be symmetric 
around x = xq. These properties can be understood from the 'one-loop dominance'. 

Among the quantities that characterize the shape of distributions, let us focus on the 
width (T, which may be defined by 



" " / 

Jo 



dx {x - xo)'^f{x) 



1 



N 



tr(A^ 



2\2 



N 



tr(^M)' 



(5.6) 



We can calculate it by the perturbative expansion around the single fuzzy sphere. Since 
the classical result is zero, it starts from the one-loop contribution, which is given by 



= 21ogiV + 0(l) 



(5.7) 



at large N (See Appendix C for a derivation). The perturbative result is independent 
of a at this order, but higher order corrections may yield terms proportional to (^)" 
{n = 1,2, • • • ). The a independence observed in our Monte Carlo data can therefore be 
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understood by the one-loop dominance at large N. In Fig. 7 we plot the variance (5.6) 
obtained from simulations with N = 8, 16, 32, 48, 64. Taking account of the 0(1) term, we 
find that the large N behavior agrees very well with the one-loop result (5.7). 

5.2 Space-time picture in the Yang-Mills phase 

Let us move on to the result for the eigenvalue distribution f{x) in the Yang- Mills phase, 
which is shown in Fig. 5 (left). If the dominant configurations in this phase were given by 
the diagonal configurations (2.4), the eigenvalues of the Casimir operator would be given 
by S/i(^M '')^ (i = 1' ■ ■ ■ )^)- Therefore the eigenvalue distribution f{x) of the Casimir 
operator (5.2) would be related to the radial density distribution p(r) of the N space-time 
points x^^^ (j = 1, • . . , TV) by 



f{x) dx = p{r) 47rr^ dr , 



(5.8) 



where x = r^. Prom (5.8) we get 



p{r) 



27rr 



(5.9) 



In reality the dominant configurations are not diagonal up to the SU(A'^) symmetry, but 
we may still define the 'radial density distribution' p{r) through (5.9). In Fig. 8 we plot 
p{r) for N = 8, 16, 32 at a = 0.0 and for = 32 at a = 0.0, • • • , 0.6. For a > 0.4 we have 
to use the zero start in order to stay in the Yang-Mills phase. The distribution at a = 
converges at large TV, and this is confirmed also for a 7^ 0. 
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Figure 8: The 'radial density distribution' p{r) is plotted for N = 8, 16, 32 at a = 0.0 (left) and 
for N = 32ata = 0.0, 0.2, 0.4, 0.6 (right). 

An interesting point to note here is that the distribution p{r) has an empty region 
< r < ro, where ro = 0.8 ~ 0.9. This may be understood as follows. If there exists an 
eigenvector of the Casimir operator Q = ^^(A^)'^ whose eigenvalue is close to zero, the 
same vector is an approximate eigenvector of Ai, A2 and A3 separately with eigenvalues 
which are also close to zero. Existence of such a vector is not allowed, however, since a 
generic configuration {Aj^; p = 1,2,3} that appears in the ensemble does not commute with 
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each other. The appearance of the empty region can thus be understood by the 'uncertainty 
principle'. Note that the scale of noncommutativity is of 0(1) since {^tiF^^) ~ 0(1) in 
the Yang-Mills phase. Therefore the size vq of the empty region should be of 0(1) as well, 
which is in agreement with our Monte Carlo results. 

Except for this empty region, the distribution p(r) decreases monotonously. This 
behavior can be understood by the one-loop effective action (2.5), which gives rise to an 
attractive potential between any pair of space-time points. The rapid fall of p{r) at large r 
is consistent with the results obtained for a = [77]. Note, however, that the precise form 
of the distribution cannot be understood solely from the perturbation theory around the 
diagonal configuration (2.4), since the higher order corrections become non-negligible as the 
space-time points come closer to each other due to the one-loop potential (2.5). For instance 
the distribution p{r) is seen to expand to larger r as a increases, but this behavior does not 
agree with the naive expectation from the one- loop effective action (2.5), where the second 
term adds to the attractive potential. We may understand it from a nonperturbative point 
of view, however, since the Chern-Simons term, which can be more negative for larger A^, 
effectively enhances configurations with a larger extent. 

6. Properties of the multi fuzzy spheres 

As we described in Section 2 the model has various multi fuzzy spheres as classical solutions 
in addition to the single fuzzy sphere. The aim of this Section is to see how they appear 
in simulations and to study their properties. 

6.1 Meta-stable states in the thermalization process 
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Figure 9: The history of the action S (left) and the eigenvalues of the Casimir operator (right) are 
plotted for = 16, a = 2.0 with the zero start. The horizontal lines represent the classical results 
for the action and the radius squared, respectively. 



Let us note first that the action for the classical solution (2.12) is given by 



24 



(6.1) 



a=l 
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which becomes minimum for the single fuzzy sphere (k = 1) (Notice the constraint (2.13)). 
Therefore at the classical level the multi fuzzy spheres are expected to appear only as 
meta-stable states. ^ 

Indeed we do find them in the process 
of thermalization. We have performed sim- 
ulations at N = 16 and a = 2.0 > 
with the zero start. Fig. 9 shows the his- 
tory of the action S and the eigenvalues of 
the Casimir operator (5.2). We have in- 
troduced a short-hand notation (ni, n2, 
• • • , Uk) for the classical solution (2.12). 
One sweep is defined as updates of all the 
elements of the matrices by the heat 
bath algorithm (See Appendix A). We see 
many plateaus before we reach the single 
fuzzy sphere after 150,000 sweeps. These 
meta-stable states can be identified with 
the multi fuzzy spheres by comparing the 
action and the eigenvalues with the clas- 
sical results (6.1) and (2.15), respectively. 

Next we start the simulation with a particular multi-fuzzy-sphere configuration, and 
consider its 'life time' defined by the number of sweeps necessary for the state to decay 
into another state. In practice the 'life time' can be extracted from the history of the 
action, which has a clear plateau as the ones in Fig. 9 (left). In Fig. 10 we show the 'life 
time' of the multi-fuzzy-sphere state (3, 5) for iV = 8 as a function of a in the log-log plot. 
Although the 'life time' thus defined is a probabilistic quantity depending on the particular 
series of random numbers, we see from Fig. 10 that it obeys a power law on the average. 
The meaning of this observation shall be discussed in the next Section. 

6.2 The 'life time' of the k coincident fuzzy spheres 
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Figure 10: The 'life time' (r) of the multi fuzzy 
spheres (3, 5) for = 8 is shown as a function of 
a in the log- log plot. The straight line represents a 
fit to the power law. 
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Figure 11: The history of the eigenvalues of the Casimir operator (5.2) is plotted for N = 16, 
a = 10.0. The initial configuration is taken to be the k coincident fuzzy spheres (6.2) with k = 2 
(left), k = A (middle) and fc = 8 (right). 



*At the quantum level the issue becomes more nontrivial; see Section 6.4. 



Prom now on we discuss the properties of the k coincident fuzzy spheres 

A,, = aL^;;^^ ® Ifc , (6.2) 

which correspond to the case rii = ■ ■ ■ = Uk = n {= of (2.12). Such a configuration is of 
particular interest since it gives rise to a NC gauge theory on the fuzzy sphere with the gauge 
group of rank k. Prom the viewpoint of perturbative calculations, the k coincident fuzzy 
spheres are easier to handle than more general multi fuzzy spheres, and we have obtained 
explicit results in the Appendices B and C. Here we perform Monte Carlo simulation with 
A'^ = 16, a = 10.0 starting from the k coincident fuzzy spheres (6.2) with k = 2,4,8. Pig. 
11 shows the history of the eigenvalues of the Casimir operator. 



22 



20 



18 



T 1 

N=16,k=2 ■ 

N=16,k=4 □ 

N=16,k=8 X 
(4/3) Log a+6.8 

N=8,k=4 A 

N=8,k=2 • 

■ (4/3) Log a+4.8 



16 



14 



12 



10 I" 



8 Lj I I I I I lJ 

4 5 6 7 8 9 10 

Log a 

Figure 12: The 'life time' (t) of k coincident fuzzy spheres with k = 2,4,8 follows a universal 
power law. The plot shows log(Tfc'^) v.s. log a for N = 8, 16. The straight lines represent the fits 
to log(T k^) = I log a + c. 

We observe that the k coincident fuzzy spheres also have a finite 'life time', and it 
becomes shorter for larger k. For all k the decay process starts with the spreading of the 
eigenvalues. We have continued the runs shown in Fig. 11 until the configurations finally 
thermalize. After visiting various multi fuzzy spheres on the way, all the simulations 
starting with different k end up with the single fuzzy sphere, the final stage looking very 
much like Fig. 9. The 'life time' r can be extracted from the plateau in the history of the 
action as in the previous Section. In Fig. 12 we plot log(r k^) against log a, which reveals 
a power law 

Tocal k-^ (6.3) 

for both A = 8 and A = 16 over a huge range of a (a = 50 ~ 20000). For a different 
algorithm, the 'life time' will be multiplied by some constant which is naively expected to 
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be independent of both a and k. If that is the case, the powers obtained here should have 
some universal meaning. 

This result is in striking contrast to the naive expectation that the decay probability 
P is given by P ~ e~'^ = exp ^— const . ^^^"^ ^ and the 'life time' is given by its inverse. We 
interpret it as a consequence of certain instability as discussed in the Appendix D. 

The power law observed for a more general multi-fuzzy-sphere state (3, 5) in the previ- 
ous Section may be understood as a consequence of the instability for shifting the center of 
the spheres relatively. This is consistent with the rather smeared distribution of the eigen- 
values observed in Fig. 9 (right) as compared with the distribution for the single fuzzy 
sphere. Such instability, however, may disappear at larger N as discussed for the coinci- 
dent case in the Appendix D. Then the decay of multi fuzzy spheres will be suppressed 
exponentially. Indeed at = 32 we observed some cases where a multi-fuzzy-sphere state 
does not decay into other states during the simulations. 

6.3 The one-loop dominance for the k coincident fuzzy spheres 

Although the k coincident fuzzy spheres have a finite 'life time', we can still measure the 
observables studied in Section 4.1 and calculate the expectation value before the decay 
actually occurs. The observables can also be calculated within the perturbation theory 
around the k coincident fuzzy spheres (6.2) neglecting the zero modes (See Appendix C). 
The leading large A'' results at the one-loop level is given by (C.9), (C.ll), (C.12) and 
(C.13). 

Fig. 13 shows the results for the observables obtained from simulations starting from 
the k = 2 coincident fuzzy spheres (6.2). As in the case of the single fuzzy sphere {k = 1), 
we find that all the observables show a discontinuity at some critical point, but this time 
the critical point is given by 

a«*^=^ ~ 2.8 . (6.4) 

Above the critical point the results agree very well with the one-loop calculation. This 
suggests that the one-loop dominance holds for k coincident fuzzy spheres as well as for 
the single fuzzy sphere. 

We can also understand the observed value (6.4) for the critical point. As we discussed 
in Section 5.1, we can put a lower bound on the critical point by assuming the one-loop 
dominance for the space-time extent (-^tr(yl^)^). For general k we find a lower boTuid 
OKcr'^ > V^k, which is satisfied by (6.4). In fact from the analysis which leads to the result 
(4.9) for A; = 1, we obtain 

4P - (^) (6.5) 
at large N. For k = 2 it gives acr'^"^ c^; \/2 (|)^''^ = 2.95, which is indeed close to (6.4). 

6.4 Comparison of the one-loop effective action 

The instability of the k coincident fuzzy spheres, which seems to be responsible for the finite 
'life time' obeying the power law, may disappear at larger A'' as discussed in the Appendix 
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Figure 13: Various observables are plotted against a for N = 32,64. Tlie initial configuration is 
taken to be the k = 2 coincident fuzzy spheres (6.2), and the measurement has been made before 
they decay into other states. 

D. Therefore by Monte Carlo simulation alone, we cannot exclude the possibility that the 
k coincident fuzzy spheres become the true vacuum of the model at some value of a in 
the large N limit. Assuming the one-loop dominance, however, we may discuss this issue 
rigorously by using the one-loop effective action, which is derived in the Appendix B. By 
comparing (B.17) and (B.23) we find that the k coincident fuzzy spheres have a smaller 
effective action than the single fuzzy sphere for 

{N 1/4 
24 logfc j 
^ _ I > = ttfe • (6.6) 

The critical value is 2.17 for k = 2, and it grows very slowly with k. 

Here the existence of the 'lower critical point' (6.5) for the k coincident fuzzy spheres 
has a dramatic consequence. Note that the 'lower critical point' dcr — 2.1v^ grows much 
faster with k than a^. As a result it turns out that < Ocr for any k > 2. This means 
that all the k{> 2) coincident fuzzy spheres are actually taken over by the Yang- Mills phase 
at a where they could have a smaller effective action than the single fuzzy sphere. Thus 
we conclude that the k coincident fuzzy spheres cannot be realized as the true vacuum in 
this model even in the large N limit. 
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7. Summciry and discussions 

In this paper we have studied nonperturbative properties of the matrix model with the 
Chern-Simons term. In particular we have demonstrated that the single fuzzy sphere 
emerges dynamically for sufficiently large a, and it describes the true vacuum of the model 
in that regime. The model may therefore serve as a nonperturbative definition of a NC 
gauge theory on the fuzzy sphere. It would be interesting to study various correlation 
functions from the field theoretical point of view and to study the UV/IR mixing effects 
arising from noncommutative geometry as has been done for the NC torus [24-26]. 

The dynamical generation of the fuzzy sphere itself is an interesting phenomenon, since 
it provides us with a concrete example in which the space-time that appears dynamically 
has lower dimensionality than the original dimensionality that the model can actually 
describe. Although the mechanism may be different, our result gives more plausibility to 
the scenario for the dynamical generation of the four-dimensional space-time in the IKKT 
matrix model [3-6]. 

When a, is smaller than a critical value, the single fuzzy sphere is no more stable, 
and the large N behavior of the system is qualitatively the same as in the pure Yang- 
Mills model (a = 0). Here the space-time looks more like a solid ball with higher density 
toward the center due to the one-loop attractive potential. It would be interesting to 
study the supersymmetric case where the one-loop attractive potential is canceled by the 
fermionic contributions at the leading order of A'^. The density distribution actually has 
an empty region in the center, which can be explained by the 'uncertainty principle', but 
this picture may depend on how one defines the density distribution when the space-time 
is noncommutative. 

The fuzzy sphere phase and the Yang- Mills phase are separated by a first order phase 
transition and we observe a strong hysteresis. One of the interesting features of the phase 
transition is that the lower critical point behaves as air ~ -^^j whereas the upper critical 

point behaves as air'* ~ 0.66 at large N . The position of the lower critical point can 
be reproduced by considering the effective potential for the scalar mode on the fuzzy 
sphere [70]. 

Various multi fuzzy spheres appear as meta-stable states at sufficiently large a, but 
we cannot exclude the possibility that they stabilize eventually in the large N limit. We 
have argued, however, that the k coincident fuzzy spheres cannot be the true vacuum even 
in the large N limit. This conclusion was obtained from the comparison of the one- loop 
effective action taking account of the existence of the 'lower critical point' for each k. 

The quantum corrections around the fuzzy-sphere solutions arc found to be dominated 
by the one-loop contribution at large N despite the absence of super symmetry. This may be 
understood by the naive power counting argument as the one in Ref . [9] . It is noteworthy 
that the 'one-loop dominance' seems to hold even for meta-stable states such as the k 
coincident fuzzy spheres. Here the measurements in simulations are made before the meta- 
stable states actually decay, while the perturbative calculations are performed omitting the 
zero modes. 
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Thanks to the meta-stabihty of the k coincident fuzzy spheres, the model may serve also 
as a toy model for the dynamical generation of gauge group in matrix models. In Ref. [78] 
it was argued that the U(fe) gauge symmetry appears in the low-energy effective theory if 
the eigenvalues of form a cluster of size k. It would be interesting to examine whether 
such a clustering is really taking place in the present model at the fully nonperturbative 
level. Prom that point of view, it would be also interesting to search for a model in which 
the gauge group of higher rank is generated in the true vacuum. 

We hope that the present work demonstrates the usefulness of Monte Carlo simulations 
in revealing various interesting dynamics of matrix models associated in particular with the 
eigenvalues of the matrices. The actual dynamics may of course depend much on the model, 
and therefore it is important to extend these studies to other models. Of particular interest 
is the impact of supersymmetry. The supersymmetric extension of the model we studied in 
this paper is known [59], but the partition function diverges [76] except for the particular 
case oi N = 2 [75] . A sensible starting point may therefore be to add the three-dimensional 
Chern-Simons term (as the one in the present model) to the D = 4 supersymmetric matrix 
model studied in Ref. [79] . Another direction is to study matrix models which have higher 
dimensional fuzzy spheres as classical solutions [80-82]. The stability issTic in this case has 
not been studied yet due to technical complications, but the application of our method is 
straightforward. Extensions to various homogeneous spaces by replacing the epsilon tensor 
in the Chern-Simons term by the structure constant of more general Lie algebra may also 
be of interest. Finally we add that studying matrix models with a cubic term may be 
of relevance to supermatrix models based on super Lie algebra [81,83] as well as to the 
quantum Hall systems [84]. 
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A. The heat bath algorithm for Monte Carlo simulations 

In this Section we comment on the algorithm used for our Monte Carlo simulation. The 
crucial point is that the Chern-Simons term is linear with respect to each of Ai, A2 and 
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74.3. This allows us to extend the heat bath algorithm used for the pure Yang- Mills model 
{a = 0) in Ref. [69] to the present model in a straightforward manner. 
Following Ref. [69] we consider the action 

-5' = y 5^ tr {Ql, - 2Q^,G^,) +2NJ2 ^^K^l) + ^^e^^^tvA^A^Ap , (A.l) 

where G^n = {A^^Ai^}. We have introduced the auxiliary fields Q^i^ (1 < /v- < < 3), 
which are N x N hermitian matrices. For fi > v we define Q^i, as Qf^i, = Qu/j,- Integrating 
out the auxiliary fields Qfj,u, we retrieve the original action (2.2). 

The auxiliary fields {Qfj,u)ij can be updated by generating a hermitian matrix with an 
appropriate Gaussian weight and shifting it by Gfxv In order to consider how to update 
Ap, we extract the part of the action (A.l) depending on Ap as 

S' = 2Nii{Sp{Apf}-Nt^{TpAp) + --- , (A.2) 

where Sp and Tp are hermitian matrices defined by 

Sp = Y^{A,f, (A.3) 

Tp = J2(^f^Qpt^ + QpuA^,) -2iaY^ ep^vA^A^ . (A.4) 

Since {Ap)ij couples only with {Ap)jk., where the index j is common, we can update ^ com- 
ponents of Ap, which do not have common indices, simultaneously. The only modification 
needed to incorporate the Chern-Simons term is the second term in (A.4). 



B. The one- loop effective action 

In this Section we formulate the perturbation theory and derive the one-loop effective 
action for the single fuzzy sphere and the A; coincident fuzzy spheres, which is discussed in 
Section 6.4. 

We decompose A^ into the classical background and the fluctuation Ap^ as 

Ap, = X^ + A^, (B.l) 

and integrate over Ap^ perturbatively. In order to remove the zero modes associated with 
the SU(A?^) invariance, we introduce the gauge fixing term and the corresponding ghost 
term 

5g.f. = -^tr[Xp,Apf , (B.2) 
5ghost = -iVtr([X^,c][^^,c]) , (B.3) 

where c and c are the ghost and anti-ghost fields, respectively. The total action reads 

'S'total = -S* + Sg,{, + Sghost ) (B.4) 
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which is given expUcitly as 



S'total = S[X] + S'kin + -S'int , (B-5) 

-N tr { ([X^, X,] - iaei^^pXp^ [i^, A,] } , (B.6) 
+^ i a €/x^p tr (i^i^ip) +iVtr(c[X^,[i^,c]]) . (B.7) 



The linear terms in X^ cancel since X^ is assumed to satisfy the classical equation of 
motion (2.3). Here and henceforth we restrict ourselves to the case where the classical 
solution X^ is proportional to a. Then, by rescaling the matrices as ^ C(Afj_, ci-^ ac, 
c ^ ac, all the terms in the total action 5totai will be proportional to a^. This means that 
the expansion parameter of the present perturbation theory is 
The effective action W is defined by 



e-^ 



= J didcdce"^*°*-i , (B.8) 



and we calculate it as a perturbative expansion W = Yl^o^J' where Wj = 0{a^^^ ■'^). 
The classical part is given by Wq = S[X], which is nothing but the action (2.2) evaluated 



at the classical solution Aj^ = X^. 



In this Section we consider the solutions of the form (2.6), for which the third term in 
(B.6) vanishes. ^ Then the kinetic term (B.6) can be written as 



S'kin = iVa^tr 



^A^{VxfA^ + ciVx?c 



(B.9) 



where we have introduced the operator 

r^M=[Xi„M], (B.IO) 
which acts on the space oi N x N matrices. The one- loop term can be obtained as 

Wi = ^Tr\og{NiVi,f} , (B.ll) 
where the symbol Tr denotes the trace of such an operator. 
B.l The single fuzzy sphere 

Let us first consider the single fuzzy sphere X^ = aL^ff^^ . The classical part is given by 

Wo = --^a^7V2(Ar2 - 1) . (B.12) 



^We will consider a more general case in Appendix D. 
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In order to calculate the one-loop term, we note that the eigenvalue problem of the operator 
(Pa)^ in the present case can be solved by the matrix analog of the spherical harmonics 
(0 < / < iV — 1, —l<m<l), which form a complete basis in the space of N x N 
matrices and have the property 

{V^^fYim = a^l{l + l) Yim . (B.13) 

Similarly to the usual spherical harmonics, they satisfy 

-^tr {YilYi^m') = 6w5mm' , (B.14) 
Ylm = (-irn-m . (B.15) 

See Ref. [59] for other properties. Prom (B.13) we find that Yqq = Iat is a zero mode, 
which corresponds to the invariance of the action (2.2) under + a^jlN- However, 

we should omit this mode when we take the trace Tr in (B.ll) since all the matrices Ji^, 
c and c are supposed to be traceless. Thus we obtain the one-loop contribution as 



Wi = -Y,{2l + l)\og[NaH{l + l)\ . (B.16) 

1=1 

At large A'" the one-loop effective action is obtained as 

Vyi_ioop^iV2 (^-^ + loga + logivj , (B.17) 

where a = as/N. 

B.2 The k coincident fuzzy spheres 

Next we consider the k coincident fuzzy spheres 

= aL^f) 1^ . (B.18) 
By setting ni = • • • = = ra (= ^) in (6.1), the classical part is obtained as 

Wo = -^{n' - 1) • (B.19) 

In order to solve the eigenvalue problem of the operator (Vx)^ in the present case, we 
consider the n x n version of the matrix spherical harmonics and introduce a k x k 
matrix e^^'^^ , whose (a, b) element is 1 and all the other elements are zero. Then as a 
complete basis oi N x N matrices, we define 

Yit'"^ ^ YL ® e("'^) , (B.20) 

which has the property 

inrY^^:^'^=a'l{l + l)Y,t''^ . (B.21) 
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Here Yq^'^^ for all the (a, b) blocks arc the zero modes. The trace part "^^=1 ^o' "^ ~ -'- 



N 

2 



should be omitted as before, but not the others. We will study the remaining (k — 1) 
zero modes in Appendix D, but here let us consider only the non-zero modes. Then the 
one-loop contribution Wi is obtained as 



n-l 



1 " 

Wi = -k'^ Y^{21 + 1) log [NaH{l + 1)] . (B.22) 

1=1 

At large N the one-loop effective action (neglecting the zero modes) reads 

W^i-ioop^iV'(-^ + loga + log^j . (B.23) 

C. The one- loop calculation of various observables 

In this Section we apply the perturbation theory discussed in Appendix B to the one- 
loop calculation of various observables which are studied by Monte Carlo simulations in 
this paper. Here we take the background to be k coincident fuzzy spheres (B.18), but the 
results for the single fuzzy sphere can be readily obtained by setting A; = 1. As in Appendix 

B. 2, we omit the zero modes for A; > 2, which will be discussed in Appendix D. 

We note that the number of loops in the relevant diagrams can be less than the order 
of ^ in the perturbative expansion since we are expanding the theory around a nontrivial 
background. At the one-loop level, the only nontrivial task is to evaluate the tadpole 
explicitly. 

C. l Propagators and the tadpole 

Using the properties (B.14), (B.15) of the matrix spherical harmonics, the propagators for 
Afj, and the ghosts are given as 

ab 1=1 m=—l 
^ 71 1 / / 1 \fn 

((^).#)«). - ^ EE E j^^, (yl^,. . 

ab 1=1 m=-l \ ' ' J 

where the symbol ( • )o refers to the expectation value calculated using the kinetic term 
Skin in (B.6) only. 

Due to the symmetries, the tadpole {A^j) can be expressed as 

(i^) = cX^ (C.3) 

with some coefHcient c. Using the identity 

tr (X^{A^)) = ctr(X^X^) 

= c^(n2-l)Ar, (C.4) 



-25- 



the coefficient c can be determined by calculating the left hand side of (C.4). 
At the leading order in we have 

^tr (x^(i^)i_loop) = (tr(X^4)tr 



- ^tr(X^A^) tT\^-ia eupaA^ApAcr 
-(tr(X^4)tr (c[X,,[i,,c]]) 



(C.5) 

Using the fact that is a linear combination of (Y^^^ ^ Ifc), we can calculate (C.5) in 
a straightforward manner. After some algebra we arrive at 



tr (x^(4)i-ioop) = -^]v^("' - 1) • 
Using (C.4) we obtain 

C.2 One-loop results for various observables 



(C.6) 
(C.7) 



Using the propagator and the tadpole obtained in the previous Section, we can evaluate 
various observables easily at the one-loop level. 

The two-point function (-^tr(yl^)^) can be evaluated as 



N 



tr(^M)' 



1— loop 



1 

N 
a" 



tr(X^X^) + 2tr (x^(i;,)i-ioop) + (tr(i^)2)o 



or 



3 "ir^ 21 + 1 



E 



At large N with fixed a = yNa, we obtain 



2\ ^C? 



1— loop 

The Chern- Simons term (M) can be evaluated as 

2i 



1 



1 



4A;2 54 



(C.8) 



(C.9) 



(M)i_ioop = [tT{XpX^Xp) + 3tr (x^X^(ip)i_ioop) 



3Af 

— ;r(fi - 1) + - 1 - 
6 ^ ' a\ 



At large with fixed a 



Na, we get 
1 



(M)l_loop 



a-^ 1 
6 fc2 (5 



(C.IO) 



(C.ll) 



The observable {jjtvF^) can be calculated in a similar manner, but we can also deduce it 
from the exact result (4.11) as 



1— loop 



(C.12) 
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Combining (C.IO) and (C.12), we get 



^2 ('^)l-loop 



a 



24 fc2 



+ 1 . 



(C.13) 



Next wc calculate the variance of the eigenvalue distribution f(x) of the Casimir 
operator studied in Section 5. The first term in (5.6) can be evaluated as follows. 



1 

iV 



2\2 



1 

1-loop N 



tr + 2tr [x^^ 



+tr ({X^^ (^-')o}) + 2tr(X,X^(i^i,)o 
+(tr(X^i^X^i^))o + (tr(i^X^i^X^))o 



16 



+ ( 1- — 



2n2 
2Z + 1 



2^/(^ + 1) 



1 



(C.14) 



The second term in (5.6) is given by the square of (C.8). The classical part of cancels 
exactly as it should. Note that since the leading 0(A^2) terms in the one-loop part also 
cancel, we have to keep track of the subleading O(logiV) terms. Using the asymptotic 
behavior 



at large n, we obtain 



(a-^)i-ioop ^ 2 log 



(C.15) 
(C.16) 



C.3 An alternative derivation 



Since tri^^ ^nd M are the operators that appear in the action 5, we can obtain their 
expectation vahics easily by using the effective action calculated for the A; coincident fuzzy 
sphere in Appendix B.2. Let us consider the action 



5(/3i,/?2,a) = ATtr ( -^[A^^A^f + 



AAA . 



(C.17) 



where we have introduced two free parameters f3i and P2, and define the corresponding 
effective action by 



Then (tr(F„,,)2) and (M) can be obtained by 



N 



4 dW 



/3i=/32=l 



(M) = 



1 dW 



(C.18) 

(C.19) 
(C.20) 



/Jl=/32 = 1 
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By rescaling the integration variables as (3^ ^ A^, we find 

W{P,,P2,a) = ^{N^ -l)logl3i + W{l,l,ap;h2) ■ (C.21) 
Using the one-loop result 

W{1, 1, «)i_ioop = 2^(n' - 1) + + 1) l°g [^«'^(^ + 1)] ' (C-22) 

which follows from (B.19) and (B.22), we can reproduce (C.IO) and (C.12). 
D. Instability of the k coincident fuzzy spheres 

In the Appendix B.2 we encountered zero modes in the one-loop perturbative expansion 
around k coincident fuzzy spheres. This corresponds to the fact that the action (2.2) does 
not change up to the second order perturbation by the deformation 

A,, = 0^4") ® Ifc + 1„ ® 5///, , (D.l) 

where 611^ (// = 1,2,3) are k x k hermitian matrices. In fact the change of the action is 
given by 

5S = Nniik [-\ [5H^,5H,f + '^iae^,p5H^5H,5H)j , (D.2) 

where the symbol tr^ refers to the trace of a A: x A" matrix. Note that SS can always be made 
negative for small 6Hp by choosing its sign appropriately when [6H^, dHy] ^ 0. Thus the k 
coincident fuzzy spheres are actually not a local minimum of the classical action, although 
this does not exclude the possibility that they stabilize due to quantum effects at finite a 
in the large N limit. 

When 5H^ {ji = 1,2,3) commute with each other, 5S in (D.2) vanishes, which means 
that such a deformation gives rise to a flat direction. In fact this is due to the fact that the 
configuration (D.l) in this case satisfies the equation of motion (2.3) for finite 6H^. Let 
us consider whether the A; coincident fuzzy spheres are stable against such a deformation 
quantum mechanically. For that purpose we calculate the one-loop effective action around 
the deformed configuration. Note first that we can diagonalize SH^j^ by applying an SU(iV) 
transformation A^ ^ UA^U^ where U is of the form U = In^V with V G SU(A;). We 
may therefore consider a classical solution 

= a (lH ®lk + ln® , (D.3) 

where 

C^ = diag(c«,--- ,cif)) (D.4) 

without loss of generality. In order to make the configuration X^j^ traceless, the k x k 
matrices should also be traceless. 

k 

trfe (C^) = 5^c(f) = 0. (D.5) 

a=l 



-28- 



The second term in (D.3) corresponds to shifting each of the k coincident fuzzy spheres, 
where the three-dimensional shift vector for the a-th sphere is given by a cj^"^ (a = 1, • • • ,k). 

Let us therefore calculate the one-loop effective action around the classical solution 
(D.3) as we did for the k coincident fuzzy spheres in Appendix B. The classical part of the 
effective action is the same as the 0^ = case (B.19). At the one-loop level, the shifted 
configuration (D.3) has the same number of zero modes as the k coincident fuzzy spheres. 



Among them, (k — 1) corresponds to changing c 



(a) 



1, • • • ,k) under the constraint 



(D.5), so we do not have to integrate over them in the present context. The remaining 
k{k — 1) modes correspond to the 'noncommutative shifts' (D.l) discussed above, but we 
simply omit them as we did in the = case. 

The kinetic term (B.6) in the present case is given by 



>S'kii 



A^a^tr 



- { {Cx + Cxf V + 2 « e^Ai^ Cx]A, + c{Cx + Cxfc 



where we introduced the operators and as 

"(LW®lfe),M' 
(1„®C^),M 



(D.6) 

(D.7) 

(D.8) 



which act on the space oi N x N matrices. The one-loop term can be given as 

Wi = ]^rr tr' log [Na^ {{Cx+Cxf 5^p + 2i e^,p } 
-rr\og\NamCx + Cxf} 



(D.9) 



where the trace tr' is taken over Lorentz indices. We calculate W\ as an expansion with 
respect to C^^ using the formulae 



2 {bi,m+iYi^l-^ + 1 



hi y^"'*^ ^ 



CxY, 



{a,b) 



2i 

lil + l)Yi 



ia,b) 



where Cp'''^ = c~p^ — c'p^ and bi^rn is given by 



bi,m = + 1) - rn{m- 1) . 

After some algebra, we get 



(D.IO) 
(D.ll) 
(D.12) 
(D.13) 
(D.14) 

(D.15) 

(D.16) 
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The first term in (D.16) is nothing but the result (B.22) for the k coincident fuzzy spheres. 
The coefficient k is given by 



Note that the first term in (D.17) grows as |logn at large n. By performing the sum in 
(D.17) numerically, we find that the coefficient k changes its sign from negative to positive 
at n = 374. 

For the cases studied in Section 6.2, where n = 2,4, 8, the k coincident fuzzy spheres are 
therefore unstable for such a shift. Note, however, that the instability for the commutative 
shift is a quantum effect unlike the instability for the non-commutative shift. Since the 'life 
time' obeys the same power law up to huge values of a, we consider that the instability 
responsible for the observed power law is dominated by the non-commutative shift. 

For a more general multi-fuzzy-sphere state (3,5) studied in Section 6.1, on the other 
hand, we don't have a counterpart of the noncommutative shift. ^ The observed power law 
therefore suggests the existence of quantum instability for the commutative shift. Thus the 
instability for the non-coincident case is expected to be much weaker than the coincident 
case. This is consistent with our observation that the power "3.8" for the former is much 
larger than | for the latter. 
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